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^ ! Abstract 

We calculated the chiral condensate and the dressed Polyakov loop in the space-time R x S 
| and R x H 3 . The chiral condensate is the order parameter for the chiral phase transition, and 

the dressed Polyakov loop is the order parameter for the deconfinement phase transition. When 
there is a current mass, critical points for the chiral and deconfinement phase transitions are 
different in the crossover region. We show that the difference is changed by the gravitational 
effect. 
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. 1 Introduction 
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■ Quantum chromodynamics (QCD) in non-perturturbative region has two important properties. 

One is breaking of chiral symmetry, another is confinement. It is believed that chiral symmetry 
breaking is restored and the confinement phase transforms into the deconfinement phase at finite 
temperature and density. The underlying QCD phase transition at temperature and density is 
constructed by these. The QCD phase transition is searched at GSI, CERN, SPS, RHIC, and LHC 
(see, eg., Ref. [TJ). However, we understand little about the explicit relation between the chiral 
^ ■ and deconfinement transitions. For example, we do not have the satisfactory effective model that 

can describe the chiral and deconfinement phase transitions. In addition, order parameters for 
deconfinement in the functional method are not simpler than one for chiral symmetry breaking [2] . 

As a theoretical approach to understand the phase structure of QCD, the lattice QCD simulation 
is a powerful method. However, there are different results for the chiral and deconfinement phase 
transitions. The Wuppertal-Budapest group showed three critical temperature for the case of Nf = 
2 + 1. For the chiral phase transition of u and d quarks, the critical temperature is J*^'^ = 151 
MeV. For s quark, critical temperature is = 175 MeV. For the deconfinement phase transition, 
the critical temperature is T% = 176 MeV [3]. Thus, two phase transitions have different critical 
temperature. On the other hand, Bielefeld group showed the result that the critical temperature 
of two phase transitions coincide at T c = 154 MeV [3]. In this way, the lattice simulation still can 
not explain the association between the chiral and deconfinement phase transitions [5]. 

In this paper, we use the Nambu-Jona-Lasinio (NJL) model to investigate the chiral and de- 
confinement phase transitions. The NJL model can describe chiral symmetry properties on QCD 
[UJ[7]. However, since the NJL model lacks confinement, the deconfinement phase transition can 
not be described. On the other hand, the Polyakov loop extended NJL (PNJL) model can describe 



the chiral and deconfinement phase transitions [81 O [TUl fTT| [12| [T3] . Although the PN JL model is 
one way to study the deconfinement phase transition, we do not use the PNJL model. Thus, we 
use the NJL model in the absence of any confinement mechanism. Instead, we calculate the dressed 
Polyakov loop p3j. Since the dressed Polyakov loop has center symmetry, this is used as the order 
parameter for the deconfinement phase transition. As shown in Ref. [15| . the dressed Polyakov 
loop without any confinement mechanism behaviors like an order parameter. (The PNJL model 
with the dressed Polyakov loop is shown in Ref |16j.) In this method, the critical points for the 
deconfinement phase transition are coincident with the chiral phase transition in the chiral limit 
m = 0. By contrast, in the case of m / 0, the critical points are not coincident in the crossover 
region |15| . This shows that there exists an intermediate region. (If the chiral and deconfinement 
phase transitions are of first order, two phase transitions should appear simultaneously.) An inves- 
tigation by the Schwinger-Dyson equation at finite temperature and density also shows that there 
is the difference between the chiral and deconfinement phase transitions [17] . 

To investigate an association between the chiral and deconfinement phase transitions from 
a different standpoint, we add a gravitational effect. For example, the separation of the chiral 
and deconfinement phase transitions in strong magnetic field is shown in Ref. |18] , Similarly, we 
investigate whether the separation of the chiral and deconfinement phase transitions are caused 
by a gravitational effect. The NJL model in curved space-time was investigated in recent year 
[HI [SB dU [221 1231 [Ml [251 126], and it was found that chiral symmetry breaking is restored. Thus, the 
chiral phase transition is affected by the curvature. Similarly, the deconfinement phase transition 
also may be affected. To investigate a gravitational effect for the chiral and deconfinement phase 
transitions, we try two cases: the positive curvature space-time R x S 3 (the Einstein universe) 
[25] . and the negative curvature space-time R x H 3 (the hyperbolic space) [26]. To calculate the 
thermodynamic potential in these cases, we use the mean field approximation. 

This paper is organized as follows. In section [21 we review the NJL model in curved space-time. 
In section [3l we introduce the dressed Polyakov loop, and calculate the (j) dependent thermodynamic 
potentials in the space time R x S 3 and R x H 3 . We show numerical results in section [31 and a 
summary is found in section [5j 



2 NJL model in curved space-time 

The SU(N C ) NJL Lagrangian with u and d quarks (the number of flavors Nf = 2) takes the form, 



where N c = 3 is the number of colors, r = (ti, T2, T3) are Pauli matrices in the flavor space, and \x 
is the quark chemical potential. We do not introduce the isospin chemical potential. 

In 4-dimensional curved space, the gamma matrices are given by the following relations [27J 



£ = q(ij^ -m + ^7o)<? + 



G 




2N C 



{l^lv} = 2 5V , {7a, 7b} = 2r/ afe , 



(2) 



with the line element 



ds 2 = g^dx^dx" = rjabe^eldx^dx". (3) 



is the tetrad and 77^ is the Minkowski metric, diag(l, — 1, — 1, — 1). Using the spin connection 
u;^ b , the spinor covariant derivative is expressed as 




By these relations, the action in curved space-time is given by 



2 



s 



G 



d^xy/^g^^ -m + H7o)q + ^tICqQ? + {qrfq) 2 ]), 



2N r 



where g = det g^ u . By applying the bosonization procedure, the action © is rewritten as 

S = J (i 4 x^5(g(i7 M V M -m-a - ij 5 (r ■ it) + fij )q - ^(c 2 + n 2 ))- 
Boson fields a and iv satisfy the equations of motion, 

G G 
a = ——qq, tt = -— qi^rq. 

We use the mean field approximation, a = (a) = const 7^ and it = (n) = 0. Using 
this approximation, one obtains the partition function at zero temperature and non-zero chemical 
potential, 



(5) 



(6) 



(7) 



VqVqVaVive 



is 



- J VqVqexp 
-Z q exp 



i / d*xV=g{q(i^V>* -M + n lQ )q - —a 1 



i I d ^ x ^-9^Q° 2 ^ 



where M = m + a. Since the integral of q, q is the Gaussian integral, log Z q is 



(8) 



log Z q = log det[»7 M V -M + nj ] 

= logdet[( 7M V^) 2 - 2i/iV - fJ 2 + M 2 ] 



(9) 



The determinant is over spinor, color, flavor, and coordinate space. We ignored a normalization 
factor. Q in the positive curvature space, RxS s (the Einstein universe) and the negative curvature 
space, R x H 3 (the hyperbolic space) was calculated in Refs [25J and [26J. 
The line element in the space-time R x S 3 is 



ds 2 = dt 2 - a 2 {d6 2 + sin 2 6dVL 2 ), 



(10) 



where the radius of the Einstein universe a is related to the scalar curvature by the relation 
R = 6/a 2 , and (KI2 is the metric on the two dimensional unit sphere. (|9|) in the space-time 
R x S 3 is given by 

i-io g z g = | f ^f^dtiiogm + uf-pii+iogm-rf-p 2 ]), (11) 
J 1=0 

where Vt = J d^Xy/^g, 1 = 0,1,2,..., & = (I + 3/2)/a and E x = 
V and the degeneracy di are expressed as 



£ 2 + M 2 . The space volume 



V = 2ir 2 a 3 , di = 2(/ + 2)(/ + l) 
The line element in the space-time R x H 3 is 



(12) 



ds 2 = dt 2 - a 2 {d6 2 + sinh 2 0dfi 2 ), 



(13) 
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where the radius of the hyperboloid a is related to the scalar curvature by the relation R = — 6/a 2 . 
The scalar curvature in the Hyperbolic space is negative. @ in the space-time R x H 3 is given by 



^logZ 9 = 6| ^ dp^p)(\og[(\(p)+^-pl) + \o g [(\(p)-^f-pl]) (14) 
where p = aVX 2 - M 2 and p(p) = {p 2 + l/4)/(2vr 2 a 3 ). 



3 (f) dependent thermodynamic potential 

To find a critical point for the chiral and deconfinement phase transitions, we consider the thermo- 
dynamic potential. The thermodynamic potential is defined as 

fi = --^io g z, (15) 

where f3 = 1/T is the inverse temperature and V is a space volume. However, this thermodynamic 
potential has no the order parameter for the deconfinement phase transition. For this reason, we 
calculate the dressed Polyakov loop. 

We introduce the dual quark condensate. The dual quark condensate is defined as[ 

2tt a A 3 /■27T 



E "=-y„ -<?./„ (16) 

The 4> dependence is caused by the U(l) valued boundary condition ip(j3,x) = e^tp^d^x), and the 
Matsubara frequency, ui n = 2irT(n + l/2)(n = 0, ±1, ±2, . . .) for fermions is replaced by u n (<p) = 
2irT(n + </>/27r). Si, which is called the dressed Polyakov loop, contains the Polyakov loop. Thus, 
S+i (or S_i) is the order parameter for deconfinement. We refer to a w as to the chiral condensate 
for simplicity. 

The log terms in ([lip and (j!4p are the same as the ordinary thermodynamic potential in the 
NJL model. Although an imaginary part is caused by the <f) dependence, we consider only a real 
part approximately. This derives a result similar to the lattice simulation [15]. Due to this, using 
the formula, 

1 sinhx 
' x 2 + (a + 2n7r) 2 2x(cosh x — cos a) ' 

we can easily obtain the <f> dependent thermodynamic potential for pip and ([14|) . ([lip with the 
U(l) valued boundary condition is expressed as 

j „ oo oo 

-— RelogZ 9 = -— E d '( lo g[^^) + ^ + ^ 2 ] + log[^(0) + (^-^) 2 ]) 

" " n=-oo 1=0 

= -£X>(ifc + E ^log[l + e- 2 ^ + ^-2e-^ + ^cos^), (18) 

1=0 e=±l 



where E\ = ^J^f + (since = m + a^, M has the <fi dependence). Similarly, (fl~4"[) with the 
U(l) valued boundary condition is expressed as 
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1 fi 00 f°° 

-— RelogZ g = -- ^ / ^(p)(log[(A(p) + / u) 2 + W 2 (0)]+ i O g [(A(/9) _ /i) 2 +w 2 ((/))]) 

/"OO T"i 

= - 12 / dpp{p)(\(p) +J2 log[l + e- 2/3 ( A+e ^ - 2e- /3 ( A+e ^ cos 0]). (19) 

6=±i 2 

Instead of using the dimension of momentum p = p/a (we will omit the Re) |26j, 

1 fi f°° 1 T 

- — ]ogZ q = ~ J dp(p 2 + ^){E p + £ -log[l + e- 2 ^ + ^-2e-^ + ^cos^), (20) 



e=±l 



where £ p = Jp 2 + M 2 . 

The real parts of (f> dependent thermodynamic potentials for the space-time Rx S 3 and R x H s 



are written as 



n? = ^ " | E <W + E | ^ + e- 2/3(£i+ ^ } - 2e-^+^) cos 0]), (21) 

«=0 e=±l 



/*0O I rp 

/ *(P 2 + -rzKEp + E ^ ^ + e -2 ^ (Ep+eM) - 2 e "^+^ cos 0]). (22) 
./o 4a £=±1 I 



4 Numerical calculation 

Since (f2"Tj) and ([52]) are divergent at large i and p, we should introduce a cutoff. We introduce the 
multiplier e~ w '/ A as the cutoff in summation over I in f)21 j) (due to the present lack of experimental 
knowledge, the cutoff in the space-time Rx S 5 can not be determined), and the momentum cutoff 
A in (|22p . Furthermore, to make dimensional quantities dimensionless, we denote the notation 
fi/A 4 — > CI, a/ A — > a,T/A — >• T, ■ ■ ■ . Then, the regularized thermodynamic potentials are written 
as 



nf = ^-£f>-<M(^+ E |log[l + e- 2 ^^)-2e-^ + ^cos^), (23) 

l ' 1=0 e=±l 



hyp _ _3_ 

> ~ 2G 



[ MP 2 + + Y ? ^gf 1 + e" 2 ^+^ - 2e- /3 (^+^ cos<£]). (24) 

* J 4a 2 ^ 2 



G, R, T, p, and cr^ in this chapter and figures represent dimensionless parameters GA 2 , R/A 2 , T/A, p/A, 
and (Jfh/A. In the chiral limit, one can obtain a critical point to calculate the global minimum of 
these. When Q(M n ^ 0) has a minimum value, chiral symmetry is broken, and when £l(M n = 0) 
has a minimum value, chiral symmetry is restored. 
The gap equations are given by 

0-^ = 0. 25 
da 6 
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In the chiral limit, when M n = a n = is obtained from this equation, chiral symmetry is restored. 
Due to this, ordinary chiral condensate a n is the order parameter for the chiral symmetry. However, 
since the chiral symmetry is broken explicitly in the case of m / 0, i.e., the critical temperature is 
defined by the peak of the chiral susceptibility, 

X = w . (26) 

Similarly, using the dressed Polyakov loop, a critical temperature for the deconfinement phase 
transition is defined by 

dT v ; 

For decision on a critical temperature, we considered order 1CP 3 for a temperature. However, for 
example, as shown in Fig. [TJ the critical temperature has an error range from the determination 
of peak. (When m = 0, the peak width is very narrow.) Taking these into account, the critical 
temperature should have a numerical error of order 10 -2 ~ 10~ 3 at least. 

We use a large value for the dimensionless current mass m. If u and d quarks are used in the flat 
space-time, the dimensionless current mass is very small (see section I4T21 or [6j[7]). However, when 
setting a small dimensionless current mass, the difference between the chiral and the deconfinement 
phase transitions is small. Due to this, we use m = 0.1 for the space-time R x S" 3 and R x H 3 . 

The summation and integration in (j23|) and (|24|) converge at curvature R < 1 in our calculation. 
The multiplier e~ ull /A is small enough and the momentum cutoff is larger than the gravitational 
scale for R < 15. In addition, the chiral phase transition in the chiral limit is observed at high 
curvature, and o§ and Si connect from low curvature (R < 1) to high curvature (R > 1) smoothly. 



4.1 The positive curvature space-time R x S 3 

In the flat space-time, the chiral phase transition is of the second order (the crossover) in the high 
temperature region in the case of m = (m 7^ 0), and is of the first order in the low temperature 
region in cases of m = and m 7^ 0. The deconfinement phase transition is the crossover in the 
high temperature region, and is of the first order in the low temperature region |15j . 

Similarly, we calculate the 4> dependent chiral condensate and the dressed Polyakov loop Si 
in the space-time R x S s . The value of the parameter G is taken as 10. We perform the numerical 
calculation for two cases, m = and m ^ 0. Figs. [2] and [3] show the <f> dependence of the o$ for 
m = and m = 0.1, respectively. o$ decreases with increasing curvature R. When compared to 
the change in temperature, the curvature simply lowers the value of The chiral condensate ov 
and the dressed Polyakov loop Si in cases of m = and m = 0.1 are shown Fig. HI Fig. H] shows 
that the chiral phase transition in the case of m = (m 7^ 0) is of the second order (crossover), 
and the deconfinement phase transition is the crossover in case of m = and m/0. On the other 
hand, since the chiral condensate and the dressed polyakov loop have discontinuous values at the 
critical points in the low temperature and high chamical potential region in the same way as the 
flat space-time, these phase transitions are of the first order. 

The critical points for the chiral phase transition are identical with the deconfinement phase 
transition in the case of m = 0, even there is the curvature (Fig. [5]). By contrary, Fig. [6] shows that 
the critical points for the chiral phase transition are not identical with the deconfinement phase 
transition in the case of m = 0.1. To clarify the difference, G is taken as 15 in Figs. [6] and [71 The 
NJL model in the ordinary four dimensional flat space-time includes this property. However, the 
difference gets smaller by increasing curvature. In addition, when m = 0.4, the difference is 0.14 at 
R = 0.1, and is 0.115 at R = 15. Thus, the difference gets larger by increasing the current mass. 
On the other hand, since the change of the difference by increasing curvature is very small and 
the peak widths of % an d t are not narrow compared to the change (see Fig. [TJ, our result bears 
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uncertainty. However, the change of the difference by increasing curvature is identified within the 
range of our numerical error. For this reason, the difference should change by increasing curvature. 

Note that the points at around <ft = ir/2 in Figs. [2] and [3] are not invariant under changes in T 
and [i. Actually, these points depend to T and fi weakly (see Fig. [8|). The contribution terms in 
<70 for 4> = 0, vr/2, ir are given by 



eP(E-n) _ i 


eP{E+y) _ 1 ' 


1 


1 


e 2P(E-n) + I 


e 2/3(E+/i) + I 


1 


1 


eplE-y) + i 


e/3{E+n) + 1 ' 



Thus, due to the factor 2, when compared to <jq and a n , may not change significantly under 
changes in T and fi. 



4.2 The negative curvature space-time R x H 3 

Chiral symmetry is broken for an arbitrary small coupling constant in the space-time R x H 3 , and 
the chiral condensate depends nonanalytical on the curvature |26j . However, we consider only the 
case of strong coupling. G is taken as 15. Figs. [9] and [10] show the <fi dependence of the for 
m = and m = 0.1, respectively. Since increases with increasing curvature \R\, the negative 
curvature yields the opposite effect. The chiral condensate a n and the dressed Polyakov loop Si in 
cases of m = and rn = 0.1 are shown Fig. [TTJ 

Also the negative curvature does not make the difference between the critical points in the case 
of m = (Fig. [T2]) . By contrary, in the case of m = 0.1, the negative curvature increases the 
difference (Figs. [131 and [bi|) . 

Unlike the space-time Rx S 3 , since the thermodynamic potential (|24p in the space-time Rx H 3 
corresponds to one in the flat space-time, one should be able to use the cutoff A in the flat space- 
time. The cutoff A in the flat space-time is estimated in the usual manner [7j. For example, 
using the pion mass 135 MeV and pion decay constant 92.4 MeV, we obtain a n A = 392 MeV with 
G = 14.6, mA = 5.6 MeV, and A = 588 MeV in the flat space-time. However, since we used the 
large value for the current mass, the current mass does not correspond to these, and is not realistic. 
Thus, we show values for the purpose of reference. Using G = 15 and m = 0.1, since we obtain 
A = 317.3 MeV in the flat space-time, at fi = 0, the difference of temperature in the space-time 
Rx H 3 is 8.2 MeV at \R\ = 0, 9.2 MeV at \R\ = 1, and 14.9 MeV at \R\ = 10. 



5 Summary 

In this paper, using the NJL model in the curved space, we investigated whether the chiral and 
deconfinement phase transitions are separated by the gravitational effect. Then, we used the dressed 
Polyakov loop as the order parameter for deconfinement. We tried the positive curvature space- 
time R x S 3 , and the negative curvature space-time R x H 3 . In both cases, the critical points are 
not different in the chiral limit m = 0. By contrast, the critical points are different for the case 
of m in the crossover region. The NJL model in the flat space has this property. However, 
the difference decreases (increases) with increasing curvature in the space- time Rx S 3 (Rx H 3 ). 
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Thus, the gravity should have a different effect on the ordinary chiral condensate and the dressed 
polyakov loop in the crossover region. 

Thus, the difference for the critical points in the flat space is caused by the presence or absence 
of the current mass m, and the gravity induces changes of the difference. In addition, a mass should 
relate to the gravity. For this reason, the ordinary chiral condensate and the dressed polyakov loop 
have a different dependence for the current mass in the crossover region (the difference increases 
with the increasing current mass in the flat space), and the positive (negative) curvature reduces 
(enhances) the contribution of the current mass. Incidentally, since the difference with increasing 
curvature for small values of m is very small, light particles are scarcely affected by the gravitational 
effect. Thus, the gravitational effect should become effective for heavy particles, such as s quark. 

Due to using the model having no a gluon, only a quark is taken into account in this paper. 
However, QCD has the gluon dynamics, and a gluon relates to confinement. Thus, to get more 
understanding of the gravitational effect for the chiral and deconfinement phase transitions, we 
must investigate QCD. In particular, the gluon may be affected by the gravity. 
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Figure 1: x and r at around peaks. The solid line indicates —0.1%, and the dotted line indicates 
t. G = 10, m = 0.1, R = 10, and /z = 0.5 
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Figure 2: (f> dependence of the in the case of m = 0. From the top to the bottom the solid lines 
correspond to R = 12, R = 13.3, and R = 14 at T = 0.5 and fi = 0.5, respectively. The dashed line 
and the dotted line correspond to T = 0.65, \i = 0.5 and T = 0.5, \i = 0.8 at R = 12, respectively. 
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Figure 3: </> dependence of the o$ in the case of m = 0.1. From the top to the bottom the solid 
lines correspond to R = 5 and R = 7 at T = 0.5 and \i = 0.5, respectively. The dashed line and 
the dotted line correspond to T = 0.8, fi = 0.5 and T = 0.5, \i = 0.8 at R = 7, respectively. 
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Figure 4: The chiral condensate a n and the dressed Polyakov loop Si in cases of m = and m = 0.1 
at R = 10 and fj, = 0.5. 




Figure 5: T — fi phase diagram for the case of m = 0. The solid line indicates the first order phase 
transition, and the dashed line indicates the second order phase transition at R = 10. 




Figure 6: T — jjl phase diagram for the case of m = 0.1. The solid line indicates the first order phase 
transition, and the dashed and dash-dotted line indicate the crossover at R = 10. The dashed line 
indicates the chiral condensate, and the dash-dotted line indicates the dressed polyakov loop. 
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Figure 7: Circle marks indicate the difference T c d — T c x . G = 15, m = 0.1, and /x = 0. 
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Figure 8: (m = 0) at around <\> = tt/2. The point at (f> = tt/2 is not invariant under changes in 
T and /i. 
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Figure 9: <fi dependence of the in the case of m = 0. From the top to the bottom the solid lines 
correspond to \R\ = 0.5 and \R\ = 0.2 at T = 0.2 and fi = 0.1, respectively (R < 0). The dashed 
line and the dotted line correspond to T = 0.25, \i = 0.1, and T = 0.2, /j, = 0.3 at \R\ = 0.5, 
respectively. 
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Figure 10: <fi dependence of the in the case of m = 0.1. From the top to the bottom the solid 
lines correspond to \R\ = 0.2 and \R\ = 0.5 at T = 0.2 and \i = 0.1, respectively (R < 0). The 
dashed line and the dotted line correspond to T = 0.25, fj, = 0.1 and T = 0.2, fj, = 0.3 at \R\ = 0.5, 
respectively. 
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Figure 11: The chiral condensate a n and the dressed Polyakov loop Si in cases of m = and 
m = 0.1 at R = -0.5 and fi = 0.1. 
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Figure 12: T — fi phase diagram for the case of m = 0. The solid line indicates the first order phase 
transition, and the dashed line indicates the second order phase transition at R = —0.5. 
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Figure 13: T — fi phase diagram for the case of m = 0.1. The solid line indicates the first order 
phase transition, and the dashed and dash-dotted lines indicate the crossover at R = —0.5. The 
dashed line indicates the chiral condensate, and the dash-dotted line indicates the dressed polyakov 
loop. 




Figure 14: Circle marks indicate the difference — T*. G = 15, m = 0.1, and \i = 0. 
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